Orientation dependence of the elastic instability on strained SiGe films 
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At low strain, SiGe films on Si substrates undergo a continuous nucleationless morphological evolution known 
as the Asaro-Tiller-Grinfeld instability. We demonstrate experimentally that this instability develops on 
Si(001) but not on Si(lll) even after long annealing. Using a continuum description of this instability, we 
determine the origin of this difference. When modeling surface diffusion in presence of wetting, elasticity and 
surface energy anisotropy, we find a retardation of the instability on Si(lll) due to a strong dependence of 
the instability onset as function of the surface stiffness. This retardation is at the origin of the inhibition of 
the instability on experimental time scales even after long annealing. 
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The complexity of electronic devices has grown exten- 
sively during the last decades. The fundamental lim- 
itations for further device improvement challenge the 
physical mechanisms at stake in the nanoscales. In 
CMOS technology, silicon-germanium (SiGe)-based ma- 
terials have shown to overwhelm silicon-only materials 
and to compete seriously with the high performances 
of III- V compounds (cut-off frequencies, current gain...). 
The understanding of the morphological evolution of thin 
SiGe films thus displays a particular interest for potential 
application. 

The coherent epitaxy of Sii-zGez films on Si sub- 
strates follows different paths depending on the growth 
parameters: 1,2 a nucleationless morphological instability 
at low x followed by island formation on top of a wetting 
layer, 2D/3D island nucleation at large x after the com- 
pletion of a thin wetting layer, step bunching on vicinal 
substrates, or plastic relaxation via nucleation of disloca- 
tions above a critical height h D . We are interested here by 
the understanding of the elastic instability onset on nom- 
inal substrates 5 and of its basic ingredients. The latter is 
reminiscent of the Asaro-Tiller-Grinfel'd (ATG) instabil- 
ity first explained in Refs. 4-6 by an enhanced surface dif- 
fusion driven by a combination of elastic relaxation and 
surface energy reduction. In addition, for thin films, wet- 
ting interactions also come into play and stabilize a 2D 
wetting film so that the instability can not occur below 
a thermodynamic critical height h c .' Moreover, it was 
shown that as-grown films on different orientations un- 
dergo different evolutions: s while the instability is man- 
ifest on Si(001), it is inhibited on Si(lll) where longer 
growth eventually leads to dislocations. We investigate 
here the influence of surface energy anisotropy on the dy- 
namics of the instability through annealing experiments 
that we confront to the continuous theory. Considering 
the extra ingredient of the instability modelling which 
is the film surface energy anisotropy, we argue that the 
difference between the two orientations is attributed to 
a retardation of the instability on Si(lll). The latter is 
characterized by a strong dependence of the instability 



onset as function of the surface stiffness. 

Growth experiments were performed by Molecular 
Beam Epitaxy (MBE) in a Riber system with a base 
pressure of ~T0-11 Torr. After 950°C in situ flashing, a 
Si buffer layer 50 nm thick is deposited at 750°C to pro- 
vide a perfectly clean reproducible flat surface. Si flux 
is obtained from an electron beam evaporator and main- 
tained constant during the growth at ~0.1ML/s. Ge is 
deposited from an effusion cell. The growth temperature 
is 550°C and the Sio.85Geo.15 composition corresponds to 
a low misfit (~0.6%) where the elastic instability occurs 
during growth without onset of island formation. Sili- 
con substrates are rotated during the experiments and 
their temperature is real time recorded. Morphologies 
z = h(x,y) are analyzed after growth by atomic force mi- 
croscopy (AFM) in air in contact and no contact modes. 
These evolutions were investigated on Si(001) and Si(lll) 
just after growth and after in situ thermal annealing at 
550°C (up to 20 h). In the present work the thickness 
(130 nm) and composition (a; = 0.15) of the SiGe layers 
were adjusted to maintain the structures below the criti- 
cal thickness of dislocation while providing enough strain 
energy to develop the instability. 
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FIG. 1. (color online) AFM top-view images of 130 run-thick 
as-grown Sio.s5Geo.15 layers (a) on (001) Si substrate and (b) 
on (111) Si substrate; side views are shaped in (c) and (d) 
respectively. 

The surface of 130 nm thick Sii-zCe^ films just after 
growth on Si(lll) and Si(100) are displayed in Fig. 1. 
On Si(001), the elastic instability is already fully devel- 
oped at the end of the lh long growth, and produces a 
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periodic corrugation characterized by a root mean square 
roughness w= ((h — {h}) 2 ) 1 ^ 2 equal to w= 1.87 nm. On 
the contrary, the experiment performed in the same ex- 
perimental conditions on a Si(lll) substrate leads to an 
almost flat surface with an intrinsic noise characterized 
by w = 0.33 nm. To investigate the influence of kinetics 
on the onset of the instability, 20 h long annealing of the 
2D film deposited on Si(lll) was performed which result 
is displayed in Fig. 2. Annealing did not allow for the de- 
velopment of the instability, since the resulting film is not 
significantly changed and is characterized by a roughness 
with u) = 0.32nm. 
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FIG. 2. (color online) (a) AFM top-view image of 130nm- 
thick Sio.85Geo.15 layers on (111) Si substrate after 20 h an- 
neal; Side view is shaped in (b). 

To give some insights into the reported difference, we 
turn to the modelling of the ATG instability which is 
driven by the elastic relaxation enforced by a surface cor- 
rugation. Its description involves a continuum represen- 
tation of surface mass currents, 4 ' where the film height 
h{x, y, t) evolves due to gradients in the chemical poten- 
tial. The latter includes a surface energy j(h,n) which 
dependences on h and on the local normal to the sur- 
face n reflect respectively wetting interactions and crys- 
talline anisotropy. :l For simplicity, we consider r y(h,n) = 
7/ [1 + 7fc(/i) + 7n(ji)] where 7/ is a reference surface en- 
ergy. The elastic chemical potential may be computed 
in the small-slope approximation in terms of the system 
Green function and of the difference between the film and 
substrate lattice parameters a^ s . 9 One can then define 
the length scale l =7/(1 - f)/[2(l + u)E(l - a 1 /a 8 ) 2 ] 
and time scale to = lo/ '(-D7/), where v and E are the film 
Poisson ratio and Young modulus, and D, the surface dif- 
fusion coefficient. The evolution equation at first order 
in h is then 

g = -A { 7(n , ft)AA + ^-|l +„„(,,)}, 

(1) 

with i,j = x,y and where the first and second term in the 
bracket describe the anisotropic surface energy, and the 
third term, the wetting interactions. The elastic energy 
density is given by the Hilbert transform T-La which acts 
in Fourier space as a multiplication by |fe|. 9 

The solution for h, with average ho, is given in Fourier 
space by h(k, t) = h(k, 0) exp(cri), with the growth rate 

a(k; n , ho) = -fe(/*o) k 2 + |fc| 3 - j(n Q ) k\ (2) 



where the first and last term describe the stabilizing 
wetting interactions and surface energy, and the second 
term, the long-range destabilizing elastic interactions. 
Anisotropy is embedded in the relative stiffness j(n) = 
[7 + d 2 j/dh 2 + d 2 "//dh 2 ] /ff, which is computed in (2) 
at no, the substrate orientation. 1 " Finally, we fit atom- 
istic calculations 11 with ^h{h) = Cw exp(— h/a* ), where, 
the variation for Ge on Si gives the extrapolation c w = 
0.05 for x = 0.15. In order to understand the dramatic 
difference of the film evolution between (001) and (111) 
orientations, we focus on the influence of the surface stiff- 
ness on the instability onset, which is known to differ 
markedly between these orientations. 

To study the thermodynamic stability of the film, we 
first consider the critical height h c above which the mor- 
phological evolution may happen as the energy of the cor- 
rugated layer becomes lower than the one of a flat layer. 
We define h* above which cr(fc) given in Eq. (2) displays 
a local maximum beside k — 0, and which depends log- 
arithmically on 7 due to the exponential dependence of 
7^. We then compute h c numerically by searching for the 
positiveness of a at this local maximum. 1 -'The value of h c 
as function of the surface stiffness is displayed in Fig. 3 
and is characterized by a logarithmic behavior with little 
variation over reasonable values of 7. Given the large de- 
posited film thickness (above 200 ML), we conclude that 
the inhibition of the elastic instability on Si(lll) cannot 
be attributed to an energetic effect linked to an enhanced 
stabilization of the dispersion relation by the surface stiff- 
ness. 
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FIG. 3. Critical height (in ML) above which the elastic insta- 
bility may develop as function of the surface stiffness 7. 

Considering this conclusion, we now study the dynam- 
ics of the instability. We look for the time t d where the 
transition between the flat fully-strained and the cor- 
rugated layers experimentally occurs. We characterize 
the emergence of the instability by the criterion that the 
roughness is greater than some value w s . Given the so- 
lution (2) for h at linear order, the roughness at time t is 
given by w 2 = (2n/L) J d 2 k\hi(k,0)\ 2 cxp(2at), where 
L is the system size in the x and y directions. We con- 
sider for the surface initial condition, a white noise of 
roughness wq distributed on L 2 modes. The criterion for 
the instability to be observable is then merely 

w 2 Jwl = (2tt)- 2 / d 2 ke 2 ^ 1 , (3) 
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which can be numerically computed, leading to the time 
t d for the onset of the instability. 




FIG. 4. Characteristic time for the instability to be fully 
developed as function of the surface stiffness. 

The resulting t d for a thick film as function of the sur- 
face stiffness is plotted in Fig. 4, considering the thresh- 
old w s /wq = 5. We find that t d increases markedly over a 
small increase of 7. becoming rapidly beyond experimen- 
tal reachable time scales for relatively modest stiffnesses, 
revealing a notable retardation of the instability dynam- 
ics induced by the surface stiffness. The inhibition of the 
elastic instability on Si(lll) may hence be attributed to 
this stiffness- induced retardation as the (111) orientation 
is known to be stiffer than the (001) orientation. 13-16 This 
higher stiffness of Si(lll) is due both to a higher step for- 
mation energy, between 0.01 eV and 0.15 eV for the two 
different step edges on Si(001) and 0.19 eV on Si(lll), 
and lower step diffusivity (inversely proportional to the 
step stiffness) 17 lA on Si(lll) and 15A on Si(001). As 
shown in Fig. 4, the evolution of t d with 7 is so large that 
the system may reach during growth the critical thickness 
for dislocation nucleation before the onset of the morpho- 
logical evolution. Hence, we attribute the transformation 
during growth on Si(lll) 8 of 2D layers into dislocated 2D 
layers to this phenomena. 

Our conjecture about the stiffness-induced retardation 
is also supported by the evolution of SiGe layers on vici- 
nal substrates. 18 On vicinal Si(001) substrates the ATG 
instability follows the same morphological evolution than 
on nominal substrate, but undergoes a morphological 
change with the transformation of square based ripples 
into ID elongated ripples. 18 On the other hand, on vic- 
inal Si(lll) (when the angle exceeds 1°), the onset of a 
step bunching instability at low misfits 19 which relaxes 
strain through edge effects is observed as opposed to the 
case on nominal substrate reported above. Considering 
the strong effect of the atomic steps on Si(lll) which en- 
force a much lower stiffness of the vicinal surface as com- 
pared to the nominal one, we advocate that the surface 
stiffness is the essential parameter ruling the dynamics of 
strained films. 

If a crucial difference between the (001) and (111) ori- 
entations concerns their stiffnesses, a relevent ingredient 
which may also come into play is the difference in their 
diffusion coefficients D. Indeed, the time scale to of the 



instability is inversely proportional to D and a kinetic 
difference could also originate from this coefficient. How- 
ever, considering the quantitative differences of the diffu- 
sion coefficients between these two orientations reported 
in Ref. 17, the inhibition of the instability on Si(lll) 
can not be explained by this effect. Indeed, the time 
constants for diffusion on Si(001) and Si(lll) at 950°C 
are 1.510 _4 s 20 and 10 _6 s JJ respectively, which testify a 
higher surface diffusion on Si(lll) compared to Si(001). 
Finally, other effects such as alloying are unlikely to have 
a strong enough influence on the instability dynamics. 

As a conclusion, we studied the orientation dependence 
of the morphological ATG instability during MBE of thin 
SiGe strained films in coherent epitaxy on Si. While the 
strain-induced corrugation is fully developed on Si(001) 
substrates, its counterpart on Si(lll) appears to be in- 
hibited even after a 20 h long annealing. We revisited 
the basic ingredients of the elastic instability and con- 
cluded that the surface stiffness has a dramatic influence 
as it enforces a significant retardation of the dynamics. 
Considering the time necessary to obtain a large enough 
roughness, we find that the instability may be postponed 
over experimentally unreachable timescales due to the 
surface stiffness and we argue that this effect causes the 
absence of the morphological evolution on Si(lll). 

The authors thank R. Kern for fruitful discussions. 
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